IAF-89-361 



Optimal Trajectories 

for the Aeroassisted Flight Experiment 


A. Miele, T. Wang, W. Y. Lee, and Z. G. Zhao 
Rice University, Houston, Tex^ 

/) 


40th Congress of the International 
Astronautical Federation 

October 7-13, 1989 
Malaga, Spain 


{'IASA-C'’-lb 6 13 3 ) 
THL At S T bD 

yn i v. > - r> n 


uPTIMAL T° A JL r T ! ’Kl l c > f'-'R 
p l T GH T fc X P P KI Mr NT ( p i <". e 

C r - C t ? 2 A 



Unc1 as 

f, 3 / 1 7 02 5 7147 



IAF-89-361 


Optimal Trajectories 

for the Aeroassisted Flight Experiment 


A. Miele, T. Wang, W. Y. Lee, and Z. G. Zhao 
Rice University, Houston, Texas 


40th Congress of the International 
Astronautical Federation 


October 7-13, 1989 
Malaga, Spain 





1 


OPTIMAL TRAJECTORIES 

12 3 

FOR THE AEROASSISTED FLIGHT EXPERIMENT ' ' 

A. Miele 4 , T. Wang, 5 , W.Y. Lee 6 , and Z. G. Zhao 7 
Rice University, Houston, Texas 


Abstract . This paper deals with the 
determination of optimal trajectories for 
the aeroassisted flight experiment (AFE). 
The intent of this experiment is to sim¬ 
ulate a GEO-to-LEO transfer, where GEO 
denotes a geosynchronous Earth orbit and 
LEO denotes a low Earth orbit. Speci¬ 
fically, the AFE spacecraft is released 
from the space shuttle and is accelerated 
by means of a solid rocket motor toward 
Earth, so as to achieve atmospheric entry 
conditions identical with those of a 
spacecraft returning from GEO. During the 
atmospheric pass, the angle of attack is 
kept constant, and the angle of bank is 
controlled in such a way that the following 
conditions are satisfied: (a) the atmos¬ 
pheric velocity depletion is such that, 
after exiting, the AFE spacecraft first 
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ascends to a specified apogee and then de¬ 
scends to a specified perigee; and (b) the 
exit orbital plane is identical with the 
entry orbital plane. The final maneuver, 
not analyzed here,includes the rendezvous 
with and the capture by the space shuttle. 

In this paper, the trajectories of an 
AFE spacecraft are analyzed in a 3D-space, 
employing the full system of 6 ODEs de¬ 
scribing the atmospheric pass. The atmos¬ 
pheric entry conditions are given, and the 
atmospheric exit conditions are adjusted in 
such a way that requirements (a) and (b) 
are met, while simultaneously minimizing 
the total characteristic velocity, hence 
the propellant consumption required for 
orbital transfer. Two possible transfers 
are considered: (IA) indirect ascent to a 
178 NM perigee via a 197 NM apogee; and 


4 

Foyt Family Professor of Aerospace 
Sciences and Mathematical Sciences, 
Aero-Astronautics Group, Rice 
University, Houston, Texas. 

5 Senior Research Scientist, Aero- 
Astronautics Group, Rice University, 
Houston, Texas. 

5 Post-Doctoral Fellow, Aero-Astronautics 
Group, Rice University, Houston, 

Texas. 

7 Research Associate, Aero-Astronautics 
Group, Rice University, Houston, 

Texas. 


ORIGINAL, PAG£ S3 

OF POOR QUALITY 



2 

(DA) direct ascent to a 178 NM apogee. 

For both transfers, two cases are 
investigated: (i) the bank angle is contin¬ 
uously variable; and (ii) the trajectory 
is divided into segments along which the 
bank angle is constant. For case (ii), the 
following subcases are studied: two seg¬ 
ments, three segments, four segments, and 
five segments; because the time duration 
of each segment is optimized, the above 
subcases involve four, six, eight, and ten 
parameters, respectively. 

It is shown that the optimal traject¬ 
ories of cases (i) and (ii) coalesce into 
a single trajectory: a two-subarc traject¬ 
ory, with the bank angle constant in each 
subarc (bang-bang control). Specifically, 
the bank angle is near 180 deg in the 
atmospheric entry phase (positive lift 
projection phase) and is near 0 deg in the 
atmospheric exit phase (negative lift pro¬ 
jection phase). It is also shown that, 
during the atmospheric pass, the peak 
values of the changes of the orbital in¬ 
clination and the longitude of the ascend¬ 
ing node are nearly zero; hence, the peak 
value of the wedge angle (angle between the 
instantaneous orbital plane and the initial 
orbital plane) is nearly zero. This means 
that the motion of the spacecraft is nearly 
planar in an inertial space. 

Key Words . Flight mechanics, hyper¬ 
velocity flight, atmospheric flight, 
optimal trajectories, aeroassisted flight 
experiment, aeroassisted orbital transfer. 


1. Introduction 

The field of aeroassisted orbital 
transfer (AOT) has received considerable 
attention in recent years. See for example 
Refs. 5-9 and references therein. In this 
paper, which is based on Refs. 1-4, the 
problem of the optimal trajectories of an 
AFE vehicle is formulated, solved, and 
analyzed. 

The aeroassisted flight experiment 
(AFE) refers to the study of the free 
flight of an autonomous spacecraft, 
shuttle-launched and shuttle-recovered. Its 
purpose is to gather atmospheric entry 
data for use in designing aeroassisted 
orbital transfer vehicles (AOTV). The 
intent of this experiment is to simulate a 
GEO-to-LEO transfer, where GEO denotes a 
geosynchronous Earth orbit and LEO denotes 
a low Earth orbit (Refs. 10-11). 

In an actual AOT transfer, GEO-to-LEO, 
the maneuver is initiated with a propulsive 
impulse at GEO so as to decelerate the 
spacecraft and force it into an elliptical 
transfer orbit leading to atmospheric 
entry. In a simulated AOT transfer, GEO- 
to-LEO, the maneuver is initiated by re¬ 
leasing the spacecraft from the space 
shuttle, which is flying at the altitude of 
160 NM above the Earth surface. By means of 
a solid rocket motor, the AFE spacecraft 
is accelerated toward Earth, so as to 
achieve atmospheric entry conditions 
identical with those of a spacecraft re¬ 
turning from GEO. Thus, for the purposes 
of this paper, the atmospheric entry 
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conditions are to be considered as given. 

During the atmospheric pass, the angle 
of attack of the AFE spacecraft is kept at 
the constant value a = 17 deg. This value 
is such that the lift L is negative. 

Control of the AFE spacecraft is achieved 
via the angle of bank \i. Hence, the pro¬ 
jected lift Lp = Lcosp (the vertical com¬ 
ponent of the lift vector) can be made 
positive or negative at will, depending on 
the value of u. 

The time history of the control, the 
bank angle p(t), is subject to two require¬ 
ments : (a) the atmospheric velocity de¬ 

pletion must be such that, after exiting, 
the AFE spacecraft first ascends to a 
specified apogee and then descends to a 
specified perigee; and (b) the exit 
orbital plane must be identical with the 
entry orbital plane; this is the same as 
stating that the exit value of the wedge 
angle must vanish; see Ref.10. Here, one 
additional requirement is considered: (c) 
the propellant consumption required for 
orbital transfer must be as small as 
possible, implying that the characteristic 
velocity AV must be as small as possible; 
in turn, minimizing AV is equivalent to 
maximizing the exit value of the horizontal 
component of the inertial velocity; this 
implies either maximizing the exit value 
of the inertial velocity or minimizing the 
exit value of the inertial path 
inclination. 

In this paper, the optimal traject¬ 
ories of the AFE spacecraft are determined 
with reference to a 3D-space and employing 


the full system of 6 ODEs describing the 
atmospheric pass. In the general formulat¬ 
ion of Refs. 1-2, the effects due to the 
rotation of the Earth and the oblateness 
of the Earth are included. In this paper, 
the former are included, but the latter 
are excluded, since they have been found 
to be small. The resulting optimal control 
problem is solved by means of the sequent¬ 
ial gradient-restoration algorithm (SGRA, 
Refs. 12-14). While this algorithm is 
available in both the primal formulation 
and the dual formulation, the former is 
used in this paper, based on previous 
experience with various AOT problems (see 
for example Refs. 15-17). 

Using SGRA, the optimal trajectories 
are computed for two possible transfers: 

(IA) indirect ascent to a 178 NM perigee 
via a 197 MM apogee; and (DA) direct 
ascent to a 178 NM apogee. For each 
transfer, two cases are investigated: (i) 
the bank angle is continuously variable; 
hence, v(t) is treated as a control; and 
(ii) the trajectory is divided into seg¬ 
ments along which the bank angle is 
constant; hence, for each segment, vi is 
treated as a parameter. For case (ii), 
the following subcases are studied: 2, 3, 

4, 5 segments; because the time duration 
of each segment is also a parameter to be 
optimized, the above subcases involve 4, 6, 
8, 10 parameters, respectively. 

For comparison purposes and only for 
Transfer (IA), a reference trajectory is 
also considered: this is a 5-segment 
trajectory, close to the nominal trajectory 
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given in Ref. 10. 

To sum up, this paper considers ten 
optimal trajectories, five for Transfer 
(IA) and five for Transfer (DA). These 
optimal trajectories are compared with 
(and are found to be superior to) the 
reference trajectory in terms of the main 
quantities of interest, namely: total 
characteristic velocity, peak heating rate, 
peak dynamic pressure, peak change of 
orbital inclination, peak change of the 
longitude of the ascending node, and peak 
value of the wedge angle. 

Content . Section 2 contains the notat¬ 
ions. The system description is given in 
Section 3, and the optimization problem is 
formulated in Section 4. The experimental 
data are given in Section 5, and the 
numerical results are presented in Section 
6. Finally, the conclusions are given in 
Section 7. 

2. Notations 

Throughout this paper, the following 
notations are employed: 

C D = drag coefficient; 

C T = lift coefficient; 

1J 

D = drag, N; 

2 

g = local acceleration of gravity, m/sec ? 
h = altitude, m; 
i = orbital inclination, rad; 

L = lift, N; 
m = mass, kg; 

r = radial distance from the center of 
the Earth, m; 

r e - radius of the Earth, m; 
r^ = radius of the outer edge of the 
atmosphere, m; 


S = reference surface area, nw 
t = time, sec; 

V = velocity, m/sec; 

V = circular velocity at r = r , m/sec; 

a a 

a = angle of attack, rad; 
y = path inclination, rad; 
n = wedge angle, rad; 

0 = longitude, rad; 

Vi = bank angle, rad; 

Vi = Earth's gravitational constant, 

3 . 2 

m /sec ; 

p = air density, kg/m'*; 

T = final time, sec; 

0 = latitude, rad; 

X = heading angle, rad; 

w = angular velocity of the Earth,rad/sec; 
Q = longitude of the ascending node, rad; 
AV - characteristic velocity, m/sec. 
Subscripts 

0 = entry into the atmosphere; 

I = exit from the atmosphere; 

II = perigee for Transfer (IA) or apogee 

for Transfer (DA); 

22 = apogee for Transfer (IA). 

Superscripts 

• = derivative with respect to time; 

= quantity defined in an inertial system 
Acronyms 

AFE = aeroassisted flight experiment; 

AFEV = aeroassisted flight experiment 
vehicle; 

AOT = aeroassisted orbital transfer; 

AOTV = aeroassisted orbital transfer 
vehicle; 

DA = direct ascent; 

2 

DP = dynamic pressure, N/m ; 

GEO = geosynchronous Earth orbit; 
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HEO = high Earth orbit; 

2 

HR = heating rate, W/m ; 

IA « indirect ascent; 

LAN = longitude of the ascending node, 

rad; 

LEO = low Earth orbit; 

ODE - ordinary differential equation; 

OT = optimal trajectory; 

RT = reference trajectory; 

SGRA = sequential gradient-restoration 
algorithm; 

WA = wedge angle, rad. 

3. System Description 

The motion of the AFE spacecraft takes 
place partly in the atmosphere and partly 
in space. For the purposes of this paper, 
the initial point is the atmospheric entry 
point; it corresponds to GEO return 
conditions and is fixed; the final point is 
located at the altitude of 178 NM and 
corresponds to circularization into LEO. 

We consider two transfer maneuvers: (IA) 
indirect ascent to a 178 NM perigee via a 
197 NM apogee; and (DA) direct ascent to a 
178 NM apogee. 

For Transfer (IA), the key points of 
the maneuver are these:point 0, atmospheric 
entry; point 1, atmospheric exit; point 22, 
apogee (h = 197 NM) of the first post- 
atmospheric transfer orbit; and point 11, 
perigee (h = 178 NM) of the second post- 
atmospheric transfer orbit. Propulsive 
impulses are applied at two points: at the 
apogee 22 in order to raise the height of 
the perigee; and at the perigee 11 in order 
to circularize the motion. 

For Transfer (DA), the key points of 


the maneuver are these: point 0, atmos¬ 
pheric entry; point 1, atmospheric exit; 
and point 11, apogee {h = 178 NM) of the 
post-atmospheric transfer orbit. A pro¬ 
pulsive impulse is applied at only one 
point: at the apogee 11 in order to circu¬ 
larize the motion. 

For the atmospheric portion (h £ h & ) 
of the trajectory of the AFE spacecraft,we 
employ an Earth-fixed system; for the space 
portion of the trajectory (h >_ h a ) , we 
employ an inertial system; here, — 
400,000 ft = 121.9 km denotes the thickness 
of the atmosphere. For h <_ h^, we compute 
the air density using the US Standard 
Atmosphere, 1976 (Ref. 18); for h h a ,we 
assume that the air density is zero. For 
both the atmospheric portion and the space 
portion of the trajectory, we neglect the 
effects due to the oblateness of the Earth; 
we assume that the gravitational field is 
central and obeys the inverse square law. 

Atmospheric Pass . With reference to 
the atmospheric portion of the trajectory 
of the AFE vehicle, the following addition¬ 
al hypotheses are employed: (a) the at¬ 
mospheric pass is made with engine shut¬ 
off; hence, the AFE spacecraft behaves as 
a particle of constant mass; (b) the angle 
of attack is constant, a = 17 deg; (c) 
under extreme hypersonic conditions, the 
dependence of the aerodynamic coefficients 
on the Mach number and the Reynolds number 
is disregarded; (d) the sideslip angle is 
zero; hence, the side force component of 
the aerodynamic force is zero. The only 
control is the angle of bank jj . 
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Differential System . With the above 
assumptions and upon using an Earth-fixed 
system, the equations of motion include the 
kinematical equations 

0 = Vcosy cosx/rcos4), (la) 

$ =-Vcosy sinx/r, (lb) 


L = (l/2)C L (o)p(h)SV 2 , (3b) 

where the air density p depends on the 
altitude h, with 

h = r - r . (3c) 

e 

Since a is constant, both the drag coef¬ 
ficient and the lift coefficient are 


r = Vsiny, (1°) 

and the dynamical equations 

V --D/m - gsiny 
2 2 

too r(sinycos $+cosysinxcos<J>sin<}>) , (Id) 

y = (L/mV) cosu+ (V/r-g/V) cosy+2wcosxcos(}) 

2 2 

+ (w r/V) (cosycos <j>-sinysinxcos$sin<{)) , 

(le) 

X = (L/mV) sinp/cosy+ (V/r) cosycosxtan<J) 

+2w (sin<j>+tanysinxcos4>) 

2 

+ (u) r/V) cosxcoscf)sin4)/cosy . (If) 

In the dynamical equations, the symbol 
u) denotes the angular velocity of the 
Earth; terms linear inware due to the 
Coriolis acceleration; terms quadratic in 
u) are due to the transport acceleration. 
Also in the dynamical equations, the 
acceleration of gravity is given by 

g = u e /r 2 , (2) 

where \i denotes the Earth's gravitational 
constant. In addition, the aerodynamic 
forces are given by 

D = (l/2)C D (a)p(h)SV 2 , 


constant. The control of the spacecraft is 
the angle of bank \i. Equations (1) must be 
integrated, subject to (2) and (3), over 
the time interval 0 <_ t <_ t. Here, the 
initial time t = 0 corresponds to atmos¬ 
pheric entry; and the final time t = t 
corresponds to atmospheric exit. 

Transformation Relations . The follow¬ 
ing transformation relations allow one to 
pass from quantities computed in an Earth- 
fixed system to quantities computed in an 
inertial system, and viceversa: 


6 — 0 + cot, 


(4a) 

<J> = , 


(4b) 

r = r. 


(4c) 

Vcosycosx = 

Vcosycosx + u>rcos4>. 

<4d) 

Vcosysiny = 

Vcosysinx t 

(4e) 

Vsiny = Vsiny. 

(4 f) 


The first three equations are linear and 
refer to the state variables appearing in 
the kinematical equations. The next three 
equations are nonlinear and refer to the 
state variables appearing in the dynamical 
equations; in spite of the nonlinearity, 
these equations can be solved explicitly 


(3a) 
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using (5), one can determine the orbital 
inclination, the longitude of the ascend¬ 
ing node, and the wedge angle. 


to yield V, y, x i- n terms of V, y, and 
viceversa. This step is omitted for the 
sake of brevity; see Refs. 1-2. 

Orbital Elements . Once the state 
variables are known in the inertial system, 
one can compute some important orbital 
quantities such as the orbital inclination 
i and the longitude of the ascending node 
ft. These quantities are supplied by the 
relations 

cosi = coscf) cosX/ (5a) 

sin{0 - ft) = coti tan$. (5b) 

Another important quantity is the wedge 
angle n f which is the angle between the 
instantaneous orbital plane and the entry 
orbital plane. This angle is given by 

cosn=sinisini 0 cos(ft-ft^)+cosicosi 0 . (5c) 

Summary . To sum up, in the Earth- 
fixed system, the equations governing the 
atmospheric pass include the differential 
system (1) and the analytical and function¬ 
al relations (2)—(3). In this formulation, 
the independent variable is the time t, 

0 < t < t. The dependent variables include 
six state variables [0(t), <J>(t), r(t), 

V(t), y(t), x (t)], one control variable 
[p(t)], and one parameter (x). If the 
initial values of the state variables are 
given and if the bank angle program is 
prescribed, the system can be integrated 
in forward time. Then, using (4), one can 
convert the state variables computed in 
the Earth-fixed system into state variables 
computed in the inertial system. Finally, 


Initial Conditions . At atmospheric 
entry, the initial values of the state 
variables 0 0 , r Q = r a , V Q , Y Q , X 0 are 

given in the inertial system. Hence, the 
initial values of the orbital elements i^, 
ft^ are known. By definition, the initial 
value of the wedge angle satisfies the 
relation Hq = 0. Using the transformation 
relations (4), the initial values of the 
state variables 0 Q , <f> Q , r Q = r Q , V Q , y Q , 

Xq become known in the Earth-fixed system. 

Final Conditions . The final time x is 
free and is to be determined indirectly as 
the time instant at which the spacecraft 
exits from the atmosphere; hence. 


(6a) 


Application of energy conservation and 
angular momentum conservation to the exit- 
to-apogee transfer orbit yields either the 
relation 


(IA) 


2 ~2 ~2 
rZ n (2V a -V^)-2r 


22 1 


~2 2~2 2 ~ 

00 r V +r V.cos y,- 
22 a a a 1 1 1 


(6b) 


or the relation 


(DA) 


2 -2 ~2 ~2 2 ~ 2 2 ~ 
rf, (2v -vf) ”2r-, r v +r vfcos y.= 0. 
11 a 1 11 a a al 1 


(6c) 


Equation (6b) applies to the indirect 
ascent case; Equation (6c) applies to the 
direct ascent case. These relations 
guarantee that, after exiting, the space¬ 
craft ascends to the specified apogee. 
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Finally, the orbital plane at atmospheric 
exit must be identical with the orbital 
plane at atmospheric entry; hence, the 
exit value of the wedge angle must vanish, 


functional to be minimized is given by 
I = AV = AV (8a) 

with 


(6 d) 


AV 11 = V a /(r a /r ll ) ' (r a /r ll )V l cos V (8b) 


Use of the transformation relations (4) in 
conjunction with (5) allows one to trans¬ 
form the final conditions (6) to include 
quantities computed in the Earth-fixed 
system. 

4. Optimal Control Problem 

Subject to the previous constraints, 
different optimal control problems can be 
formulated, depending on the performance 
index chosen and the type of transfer 
maneuver considered. Here, we focus 
attention on the minimization of the total 
characteristic velocity AV, which is a 
measure of the propellant consumption 
required for orbital transfer. 

Problem (IA ).This problem refers to 
the indirect ascent to a 178 NM perigee 
via a 197 NM apogee. The functional to be 
minimized is given by 


I - AV = AV 22 + AV i;l , 


(7a) 


with 


AV 22= V [2r a r ll /(r ll r 22 + r 22 ): 


(r a /r 22 )V l COSY l' 


(7b) 


AV 11 = V a /[2r a r 22 /(r U r 22 + r ll> 1 


V< r a /r ll>- 


(7c) 


Problem (DA ). This problem refers to 
the direct ascent to a 178 NM apogee. The 


Alternative Formulations . Inspection 
of the performance indexes (7) and (8) 
shows that they include a constant part and 
a part which is linear in the exit value of 
the horizontal component of the inertial 
velocity = Vcosy. Hence, the minimi¬ 
zation of AV is equivalent to maximizing 
the functional 

J “ ^hl = V-jCOsy^ (9a) 

In turn, because of the final conditions 
(6b) or (6c), the maximization of the 
functional (9a) implies either maximizing 
the functional 

K = V (9b) 

or minimizing the functional 

H = y^ • (9 C ) 

Continuous vs Discrete Control . For 
both Problems (IA) and (DA), two cases are 
investigated: continuous control and dis¬ 
crete control. 

(i) Continuous Control. Here, the 

bank angle is continuously variable. Hence, 
y (t) is treated as a control variable. 

(ii) Discrete Control. Here, the 

trajectory is divided into s segments 
along which the bank angle is constant; 
hence, for each segment, the bank angle 

y^ is treated as a parameter. Because the 



time duration of each segment is also a 

parameter to be optimized, the number of 

parameters is twice the number of segments. 

The following subcases are studied: s = 2, 

3, 4, 5. These subcases involve n =4, 

P 

6, 8, 10 parameters, respectively. 

Algorithm . For both cases (i) and 
(ii)/ Problems (IA) and (DA) are optimal 
control problems of the Bolza type. They 
can be solved using the sequential 
gradient-restoration algorithm (SGRA, Refs. 
12-14) employed in conjunction with the 
primal formulation (PSGRA, Ref. 14). 

The sequential gradient-restoration 
algorithm is an iterative technique which 
involves a sequence of cycles, each cycle 
including two phases: the gradient phase 
and the restoration phase. In the gradient 
phase, the value of the augmented function¬ 
al is decreased, while avoiding excessive 
constraint violation. In the restoration 
phase, the value of the constraint error is 
decreased, while avoiding excessive change 
in the value of the functional. In a 
complete gradient-restoration cycle, the 
value of the functional is decreased, while 
the constraints are satisfied to a pre¬ 
selected degree of accuracy. Thus, a 
succession of suboptimal solutions is 
generated, each new solution being an 
improvement over the previous one from the 
point of view of the value of the function¬ 
al being minimized. 

The convergence conditions are re¬ 
presented by the relations 

p £ e l' Q £ e 2 * 


Here, P is the norm squared of the error 
in the constraints; Q is the norm squared 
of the error in the optimality conditions; 
and e^, are preselected small, positive 
numbers. 

5. Experimental Data 

The following data were used in the 

numerical experiment (for more significant 

digits, see Ref. 4). 

Earth's Physical Constants . The radius 

of the Earth is r = 6378 km; the radius 
e 

of the outer edge of the atmosphere is 
r = 6500 km; the thickness of the atmos- 

a 

phere is h = 121.9 km; the Earth's 
a 

gravitational constant is = 0.3986E+15 
3 2 

m /sec ; the circular velocity at r = r 

<3 

is V a = 7.831 km/sec; the angular velocity 
of the Earth is m = 0.7292E-04 rad/sec. 

Spacecraft Data . The mass of the AFE 

spacecraft is m = 1678 kg; the reference 

2 

surface area is S = 14.31 m ; the space¬ 
craft flies at constant angle of attack, 
a = 17.00 deg; the lift coefficient is 
C L =-0.3707; and the drag coefficient is 
C D = 1.315. 

Atmospheric Model . The assumed 
atmospheric model is that of the US 
Standard Atmosphere, 1976 (Ref. 18). In 
this model, the values of the density are 
tabulated at discrete altitudes. For 
intermediate altitudes, the density is 
computed by assuming an exponential fit 
for the function p(h) . 

Heating Rate . The stagnation point 
heating rate is assumed to be given by the 


( 10 ) 


formula 
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107 

HR = C/(p/p R ) (V/V R ) J * U/ . (11) 

Here, p R is a reference density, the density 
at 1^= 60 km; V R = V & is a reference 
velocity; the constant C represents the 
stagnation point heating rate at p = p R 

2 

and V = V and is given by C = 282.3 W/cm ; 

K 

note that the value of this constant is 
based on a nose radius of one foot. 

Atmospheric Entry Conditions . In the 
inertial system, the given initial 
conditions are as follows: the longitude 
is 0 q = -134.52 deg; the latitude is = 
-4.49 deg; the altitude is h^ = 121.9 km; 
the radius is r^ = 6500 km; the velocity 
is Vq = 10.31 km/sec; the path inclination 
is Yq = -4.49 deg; the heading angle is 
Xq = -28.13 deg; the orbital inclination 
is ig = 28.45 deg; the longitude of the 
ascending node is ftg = -126.19 deg; the 
wedge angle is n q = 0.00 deg. 

Atmospheric Exit Conditions . In the 
inertial system, the desired final 
conditions are as follows: the altitude is 
h^ - 121.9 km; the radius is r^ = 6500 km; 
the orbital inclination is i^ = 28.45 deg; 
the longitude of the ascending node is 
= -126.19 deg; the wedge angle is = 
0.00 deg. 

Transfer (IA) ♦ This is the indirect 
ascent to a 178 NM perigee via a 197 NM 
apogee. The required apogee conditions are 
as follows: the altitude is ^22 ~ 364.8 km; 
the radius is = 6743 km; and the path 

inclination is y 2 2 = 0*0° deg. The required 
perigee conditions are as follows: the 
altitude is h^ = 329.7 km; the radius is 


r^ = 6708 km; and the path inclination 
is " 0.00 deg. 

Transfer (DA ). This is the direct 
ascent to a 178 NM apogee. The required 
apogee conditions are as follows: the 
altitude is h^ = 329.7 km; the radius is 
r ll ~ 6708 km; and the path inclination 
is Yj^ = 0*00 deg. 

6. Numerical Results 

The optimal control problem formulated 
in Sections 3-4 was solved with the 
sequential gradient-restoration algorithm 
for the experimental data outlined in 
Section 5. Both Transfer (IA) and Transfer 
(DA) were optimized from the point of view 
of the characteristic velocity AV. For 
each transfer, five optimal trajectories 
(OT) were computed: the continuous control 
OT; the discrete control OT with s = 2 
segments; the discrete control OT with s = 

3 segments; the discrete control OT with 
s = 4 segments; and the discrete control 
OT with s = 5 segments. For Transfer (IA), 
a reference trajectory (RT) was also 
computed; this is a nonoptimal discrete- 
control trajectory with s = 5 segments, 
obtained by means of the restoration 
algorithm, starting from a nominal traject¬ 
ory supplied in Ref.10. Thus, a total of 11 
trajectories were computed, ten optimal 
trajectories and one reference trajectory. 
The numerical results are shown in Tables 
1-2 and Figs. 1-2. 

Table 1 shows the values of the 
characteristic velocity AV and the flight 
time t for all the trajectories computed. 
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For Transfer (IA), the reference trajectory 
requires AV = 98.0 m/sec and a flight time 
t = 488 sec; also for Transfer (IA), the 
optimal trajectories require AV = 81.9 
m/sec and a flight time t = 796 sec. For 
Transfer (DA), the optimal trajectories 
require AV = 72.0 m/sec and a flight time 
t = 820 sec. To sum up, the optimal 
trajectories are considerably more 
efficient propulsively than the reference 
trajectory and they are characterized by a 
longer flight time. 

An interesting result of the analysis 
is that, for each transfer, the optimal 
continuous-control trajectory and the 
optimal discrete-control trajectories 
coalesce into a single trajectory: a two- 
subarc trajectory, with the bank angle 
constant in each subarc (bang-bang control). 
Specifically, the bank angle is near 180 
deg in the atmospheric entry phase 
(positive lift projection phase) and is 
near 0 deg in the atmospheric exit phase 
(negative lift projection phase). This is 
true for both Transfer (IA) and Transfer 
(DA). See Figs. 1-2. 

In the light of the above statement, 
we now focus attention on only three 
trajectories: for Transfer (IA), the re¬ 
ference trajectory and the two-subarc 
optimal trajectory; for Transfer (DA), the 
two-subarc optimal trajectory. For these 
trajectories. Table 2 shows the values of 
the following quantities: the character¬ 
istic velocity; the flight time; the 
minimum altitude; the peak dynamic 


pressure; the peak heating rate; the peak 
change of the orbital inclination; the peak 
change of the longitude of the ascending 
node; and the peak value of the wedge 
angle. 

Table 2 illustrates clearly the 
advantages of the OTs with respect to the 
reference trajectory. Concerning the 
longitudinal motion, by comparison with 
the RT, the OTs are characterized by 
smaller AV and also by smaller values of 
the peak dynamic pressure and the peak 
heating rate; this is due to the fact that 
the minimum altitude of the OTs is about 3 
km higher than the minimum altitude of the 
RT, Concerning the lateral motion, by 
comparison with the RT, the OTs exhibit 
smaller values of the peak change of the 
orbital inclination and the peak change of 
the longitude of the ascending node; hence, 
the peak wedge angle of the OTs is smaller 
than the corresponding quantity for the RT. 
Indeed, it is suprising that max|Ai| and 
max|Aft| are nearly zero for the OTs; hence, 
max(n) is nearly zero for the OTs. This 
means that, for efficient flight, the 
motion of the spacecraft must be nearly 
planar in an inertial space; in other 
words, one must avoid energy dissipation 
associated with the lateral motion. 

More details on the trajectories 
computed can be found in Figs. 1 and 2. 
Figure 1 refers to Transfer (IA), and Fig. 

2 refers to Transfer (DA). These figures 
include the time histories of the following 
quantities: the bank angle (Figs. IA and 
2A); the altitude (Figs. IB and 2B); the 
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relative velocity {Figs. 1C and 2C); the 
relative path inclination (Figs. ID and 
2D); the dynamic pressure (Figs. IE and 
2E)? the heating rate (Figs. IF and 2F); 
and the wedge angle (Figs. 1G and 2G). 

Remark . With reference to the two- 
subarc OTs of Transfer (IA) and Transfer 
(DA), we observe that the value of the 
bank angle is \x = 176.7 deg for the at¬ 
mospheric entry phase and \i - 5.5 deg for 
the atmospheric exit phase. Let L denote 
the lift, and let 


L 

P 


Lcosy, 


L = Lsiny 

q 


( 12 ) 


denote the vertical component and the 
lateral component of the lift. Let C , C , 

Lt i-rp 

C T denote the associated lift coefficients, 
Lq 

and observe that 


C L P = C L COSlJ ' C Lq = C L sinU - (13) 

We recall that the AFE spacecraft flies at 
constant angle of attack a = 17 deg, 
corresponding to = -0.3707 and C D = 

1.315. With this understanding, we see that, 
for the atmospheric entry phase, 

C T = -0.998 C T , C T = 0.058 C T ; (14) 

Lp Ju LiQ Lt 

and we see that, for the atmospheric exit 
phase, 

C T * 0.995 C T , C T = 0.096 C T . (15) 

Lp L Lq L 

Therefore, the following comments arise: 

(i) in the atmospheric entry phase, 
the vertical component of the lift is 
directed upward and its modulus is only 
2/1000 smaller than the lift modulus; this 
causes the path inclination to increase 


gradually from the entry negative value to 
nearly zero value; 

(ii) in the atmospheric exit phase, 
the vertical component of the lift is 
directed downward and its modulus is only 
5/1000 smaller than the lift modulus; this 
offsets the centrifugal force effects due 
to the curvature of the Earth, so as to 
ensure exit conditions compatible with the 
desired apogee requirement; 

(iii) the lateral component of the 
lift during the atmospheric entry phase 
and the lateral component of the lift 
during the atmospheric exit phase have the 
same sign and the same order of magnitude; 
their modulus is 5.8% of the lift modulus 
in the atmospheric entry phase and 9.6% of 
the lift modulus in the atmospheric exit 
phase; these lateral lift components are 
directed in such a way that they nearly 
offset the effects due to the Earth's rot¬ 
ation, so as to ensure exit conditions 
compatible with the desired wedge angle 
requirement. 

7. Conclusions 

In this paper, the trajectories of an 
AFE spacecraft are analyzed in a 3D-space, 
employing the full system of 6 ODEs des¬ 
cribing the atmospheric pass. It is assumed 
that the angle of attack is constant and 
that the spacecraft is controlled via 
only the angle of bank. It is also assumed 
that the atmospheric entry conditions are 
given and that the atmospheric exit 
conditions are adjusted in such a way that 
the following conditions are satisfied: (a) 
the atmospheric velocity depletion is such 
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that, after exiting, the AFE spacecraft 
first ascends to a specified apogee and 
then descends to a specified perigee; and 
(b) the exit orbital plane is identical 
with the entry orbital plane. 

Under the above conditions, optimal 
trajectories are determined by minimizing 
the total characteristic velocity, hence 
the propellant consumption required for 
orbital transfer. The resulting optimi¬ 
zation problems are solved by means of the 
sequential gradient-restoration algorithm. 

Two possible transfers are consider¬ 
ed: (IA) indirect ascent to a 178 NM 
perigee via a 197 NM apogee; and (DA) 
direct ascent to a 178 NM apogee. For both 
transfers, two cases are investigated: (i) 
the bank angle is continuously variable; 
and (ii) the trajectory is divided into 
segments along which the bank angle is 
constant. For cases (ii), the following 
subcases are studied: two segments, three 
segments, four segments, and five seg¬ 
ments; because the time duration of each 
segment is optimized, the above subcases 
involve four, six, eight, and ten para¬ 
meters, respectively. For comparison 
purposes and only for Transfer (IA), a 
reference trajectory is also considered: 
this is a five-segment trajectory, close 
to the nominal trajectory given in Ref. 

10 . 

From the extensive numerical comput¬ 
ations, the following conclusions arise: 

(i) The optimal trajectories are 
two-subarc trajectories, with the bank 
angle constant in each subarc; hence, the 


control is bang-bang. 

(ii) In the atmospheric entry phase, 
the bank angle is near 180 deg, yielding a 
positive projected lift, which in turn 
causes the path inclination to increase 
gradually from the entry negative value to 
nearly zero value. 

(iii) In the atmospheric exit phase, 
the bank angle is near 0 deg, yielding a 
negative projected lift, which offsets the 
centrifugal force effects due to the 
curvature of the Earth, so as to ensure 
exit conditions compatible with the desired 
apogee requirement. 

(iv) The lateral component of the 
lift during the atmospheric entry phase 
and the lateral component of the lift dur¬ 
ing the atmospheric exit phase have the 
same sign and the same order of magnitude; 
they are directed in such a way that they 
nearly offset the effects due to the 
Earth's rotation. In this way, the instant¬ 
aneous orbital plane is almost identical 
with the initial orbital plane, meaning 
that the wedge angle is nearly zero during 
the atmospheric pass. This means that, for 
efficient flight, the motion of the AFE 
spacecraft is nearly planar in an inertial 
space; in other words, one must avoid 
energy dissipation associated with the 
lateral motion. 

(v) Comparison of the optimal 
trajectories and the reference trajectory 
shows that the OTs are superior to the RT 

in terms of the main quantities of interest, 
namely, the characteristic velocity, the 
peak dynamic pressure, the peak heating 
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rate, and the peak wedge angle. In part- 6. 

icular, for Transfer (IA), the charact¬ 
eristic velocity is AV = 98.0 m/sec for 
the RT and AV = 81.9 m/sec for the OT. For 
Transfer (DA) , the characteristic velocity 
of the OT is 72.0 m/sec. 

References 

1. MIELE, A., ZHAO, Z. G., and LEE, 

W. Y., "Optimal Trajectories for the 
Aeroassisted Flight Experiment, Part 

1, Equations of Motion in an Earth- 
Fixed System", Rice University, Aero- 
Astronautics Report No. 238, 1989. 

2. MIELE, A., ZHAO, Z. G., and LEE, 

W. Y., "Optimal Trajectories for the 
Aeroassisted Flight Experiment, Part 

2, Equations of Motion in an Inertial 
System", Rice University, Aero- 
Astronautics Report No. 239, 1989. 

3. MIELE, A., WANG, T., LEE, W. Y., and 
ZHAO, Z. G., "Optimal Trajectories 
for the Aeroassisted Flight Experi¬ 
ment, Part 3, Formulation, Results, 
and Analysis", Rice University, Aero- 
Astronautics Report No. 242, 1989. 

4. MIELE, A., WANG, T., LEE, W. Y., WANG 
H., and WU, G. D. , "Optimal Traject¬ 
ories for the Aeroassisted Flight 
Experiment, Part 4, Data, Tables, 
and Graphs", Rice University, Aero- 
Astronautics Report No. 243, 1989. 

5. WALBERG, G. D., "A Survey of Aero¬ 
assisted Orbital Transfer", Journal 
of Spacecraft and Rockets, Vol. 22, 

No. 1, pp. 3-18, 1985. 


MEASE, K. D., and VINH, N. X., 
"Minimum-Fuel Aeroassisted Coplanar 
Orbital Transfer Using Lift Modulat¬ 
ion", Journal of Guidance, Control, 
and Dynamics, Vol. 8, No. 1, pp. 134- 
141, 1985. 

7. VINH, N. X., and HANSON, J. M. , 
"Optimal Aeroassisted Return from High 
Earth Orbit with Plane Change", Acta 
Astronautica, Vol. 12, No. 1, pp. 11- 
25, 1985. 

8. MIELE, A., Editor, "Special Issue on 
Hypervelocity Flight", Journal of the 
Astronautical Sciences, Vol. 36, Nos. 
1-2, pp. 1-197, 1988. 

9. MEASE, K. D., "Optimization of Aero¬ 
assisted Orbital Transfer: Current 
Status", Journal of the Astronautical 
Sciences, Vol. 36, Nos. 1-2, pp. 7-34, 
1988. 

10. ANONYMOUS, N. N., "Aeroassisted Flight 
Experiment: Preliminary Design Docu¬ 
ment" ,NASA Marshall Space Flight 
Center, 1986. 

11. GAMBLE, J. D., CERIMELE, C. J., MOORE, 
T. E., and HIGGINS, J"Atmospheric 
Guidance Concepts for an Aeroassist 
Flight Experiment", Journal of the 
Astronautical Sciences, Vol. 36, Nos. 
1-2, pp. 45-71, 1988. 

12. MIELE, A., and WANG, T., "Primal-Dual 
Properties of Sequential Gradient- 
Restoration Algorithms for Optimal 
Control Problems, Part 1, Basic 
Problem", Integral Methods in Science 
and Engineering, Edited by 



15 


F. R. Payne et al, Hemisphere 
Publishing Corporation, Washington, 

DC, pp. 577-607, 1986. 

13. MIELE, A., and WANG, T., "Primal-Dual 
Properties of Sequential Gradient- 
Restoration Algorithms for Optimal 
Control Problems, Part 2, General 
Problem", Journal of Mathematical 
Analysis and Applications, Vol. 119, 
Nos. 1-2, pp. 21-54, 1986. 

14. MIELE, A., WANG, T., and BASAPUR, 

V. K., "Primal and Dual Formulations 
of Sequential Gradient-Restoration 
Algorithms for Trajectory Optimizat¬ 
ion Problems", Acta Astronautica, 

Vol. 13, No. 8, pp. 491-505, 1986. 

15. MIELE, A., BASAPUR, V. K., and MEASE, 
K. D., "Nearly-Grazing Optimal 
Trajectories for Aeroassisted Orbital 
Transfer", Journal of the Astro- 
nautical Sciences, Vol. 34, No. 1, 
pp. 3-18, 1986. 

16. MIELE, A., BASAPUR, V. K., and LEE, 

W. Y., "Optimal Trajectories for 
Aeroassisted, Coplanar Orbital 
Transfer", Journal of Optimization 
Theory and Applications, Vol. 52, 

No. 1, pp. 1-24, 1987. 

17. MIELE, A., BASAPUR, V. K., and LEE, 

W. Y., "Optimal Trajectories for 
Aeroassisted, Noncoplanar Orbital 
Transfer", Acta Astronautica, Vol. 15, 
Nos. 6-7, pp. 399-411, 1987. 

18. NOAA, NASA, and USAF, "US Standard 
Atmosphere, 1976", US Government 
Printing Office, Washington, DC, 1976. 


Supplementary References 

19. LONDON, H. S., "Change of Satellite 
Orbit Plane by Aerodynamic Maneuvering", 
Journal of the Aerospace Sciences, 

Vol. 29, No. 3, pp. 323-332, 1962. 

20. HULL, D. G., GILTNER, J. M., SPEYER, 

J. L., and MAPAR, J., "Minimum-Energy- 
Loss Guidance for Aeroassisted Orbital 
Plane Change", Journal of Guidance, 
Control, and Dynamics, Vol. 8, No. 4, 
pp. 487-493, 1985. 

21. CALISE, A. J., "Singular Perturbation 
Analysis of the Atmospheric Orbital Plane 
Change Problem", Journal of the 
Astronautical Sciences, Vol. 36, Nos. 

1-2, pp. 35-43, 1988. 

22. HULL, D. G., McCLENDON, J. R., and 
SPEYER, J. L., "Aeroassisted Orbital 
Plane Change Using an Elliptic Drag 
Polar", Journal of the Astronautical 
Sciences, Vol. 36, Nos. 1-2, pp. 73- 
87, 1988. 

23. MEASE, K. D., LEE, J. Y., and VINH, 

N. X., "Orbital Changes during Hyper¬ 
velocity Aerocruise", Journal of the 
Astronautical Sciences, Vol. 36, Nos. 

1-2,pp. 103-137, 1988. 

24. VINH, N. X., and LU, P., "Chebyshev 
Minimax Problems for Skip Traject¬ 
ories", Journal of the Astronautical 
Sciences, Vol. 36, Nos. 1-2, pp. 179- 
197, 1988. 


16 


Table 1. Characteristic velocity and flight time. 


Traj ectory 

Control 

s 

AV(m/sec) 

t (sec) 

IART 


Discrete 

5 

97.97 

487.6 

IAOT 


Continuous 

00 

81.86 

795.3 

IAOT 


Discrete 

2 

81.86 

795.7 

IAOT 


Discrete 

3 

81.86 

795.7 

IAOT 


Discrete 

4 

81.86 

795.8 

IAOT 


Discrete 

5 

81.84 

796.2 

DAOT 


Continuous 

OO 

72.05 

817.6 

DAOT 


Discrete 

2 

72.03 

819.5 

DAOT 


Discrete 

3 

72.03 

819.6 

DAOT 


Discrete 

4 

72.03 

819.6 

DAOT 


Discrete 

5 

72.01 

819.9 

IA = 

OT = 

s = 

indirect ascent, DA = 

optimal trajectory, RT = 

number of segments. 

direct ascent, 

reference trajectory, 



Table 2. Comparison of AFE trajectories. 


Quantity 

IART 

s = 5 

IAOT 

s = 2 

DAOT 

s = 2 

Units 

AV 

97.97 

81.86 

72.03 

m/sec 

T 

487.6 

795.7 

819.5 

sec 

min(h) 

74.6 

77.5 

77.5 

km 

max(DP) 

1629 

1174 

1174 

N/m 2 

max (HR) 

158 

148 

148 

W/cm 2 

max|Ai| 

1.08 

0.01 

0.01 

deg 

max|AH| 

0.32 

0.00 

0.00 

deg 

max(n) 

1.08 

0.01 

0.01 

deg 

IA = indirect 

ascent, 

DA = 

direct ascent. 


OT = optimal 

trajectory. 

RT = 

reference trajectory, 

s = number of segments. 








FIG.lfl. COMPRRI SON OF RFE TRRJECTORIES , 
TRANSFER (IR) . BANK ANGLE- 



FIG-2R. COMPARISON OF AFE TRAJECTORIES 
TRANSFER (DA ) . BANK ANGLE - 
















FIG.ZB. COMPARISON OF AFE TRAJECTORIES 
TRANSFER (DA), ALTITUDE. 
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FIG.1C. COMPARISON OF AFE TRAJECTORIES, 

TRANSFER (IA), VELOCITY, RELATIVE- 



FIG.2C. COMPARISON OF AFE TRAJECTORIES. 

TRANSFER (DA), VELOCITY, RELATIVE. 
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FIG.ID. COMPARISON OF AFE 
TRANSFER ( IA) . 
PATH INCLINATION. 



FIG.2D. COMPARISON OF AFE 
TRANSFER (DA). 
PATH INCLINATION. 
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